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We review recent exact results for the free expansion of impenetrable bosons on one- 
dimensional lattices, after switching off a confining potential. When the system is initially 
in a superfiuid state, far from the regime in which the Mott-insulator appears in the mid- 
dle of the trap, the momentum distribution of the expanding bosons rapidly approaches 
the momentum distribution of noninteracting fermions. Remarkably, no loss in coherence 
is observed in the system as reflected by a large occupation of the lowest eigenstate of 
the one-particle density matrix. In the opposite limit, when the initial system is a pure 
Mott insulator with one particle per lattice site, the expansion leads to the emergence of 
quasicondensates at finite momentum. In this case, one-particle correlations like the ones 
shown to be universal in the equilibrium case develop in the system. We show that the 
out-of-equilibrium behavior of the Shannon information entropy in momentum space, 
and its contrast with the one of noninteracting fermions, allows to differentiate the two 
different regimes of interest. It also helps in understanding the crossover between them. 

Keywords: nonequilibrium; hardcore; fermionization. 



1. Introduction 

Physical systems displaying the effects of a reduced dimensionality have been in- 
creasingly attracting the interest of experimental and theoretical physicists in the 
past decades. Traditionally a condensed-matter subject, it has received a lot of 
attention within the last years in the framework of ultracold quantum ga ses. T he 
realization of Bose-Einstein co ndensates (BE C) in very anisotropic traps^^^ or 
their loading on optical lattices * * * * * allow experimentalists to obtain a rich 
variety of (quasi-)one-dimensional (ID) systems where the reduced dimensionality 
rules the physics. 

In particular, the ID regime where bosons behave as impenetrable particlespEl 
also known as hardcore bosons (HCB's), has generated much interest recently. This 
regime, which has been called the Tonks-Girardeau (TG) regime,^^ can be ob- 
tained for large positive three-dimensional scattering lengths, low densities, and 



2, 2008 16:2 WSPC/INSTRUCTION FILE Expansionl 



low temperatures.E^EJ Contrary to weakly interacting bosons, HCB's in ID 
share niany properties with noninteracting spinless fermions to which they can be 
mappedM-'J (The TG regime is sometimes referred as the fermionized regime.) Prop- 
erties like the energy, density profiles, and density-density correlations, are identical 
in both systems. On the contrary, one-particle correlations, and related quantities 
like the momentum distribution function (n^), and the so-called natural orbitals, 
are very different. 

The ID TG gas, introduced by Girardeau in the 60's,^2 has been ex- 
tensively stu died ip the liter ature by means of differ ent techniqu es, both in 
homogeneous ,1^ -*- 6 | 17 | 18 | 19 | 20 | g^j^j harmonically trappeci^HlEniSl systems. In 
ID, even at zero temperature, HCB's do not exhibit BEC. As shown by Lenard,^! 
and Vaidya and Tracy^l one-particle correlations decay with a power law ~ 
1/y/x. Hence, the highest occupied effective single-particle state (the lowest nat- 
ural orbital) in the system has an occupation Aq ^ \/iVf, {N}, being the number 
of bosons), i.e., Aq ^ oo for N}, oo, but Xo/Nb — > 0, and we call this state a 
quasicondensate. A similar result has been obtained in harmonic traps. 

Once the ID regime has been achi eved exp erimentallv. for example by means of 
a two-dimensional (2D) optical latticej I I I Hhe addition of a new lattice along the 
ID axes of the tubes- facilitates the achievement of the TG regime with respect 
to the continuum case. This is because it allows experimentalists to change the 
effective mass of the particles, and consequently the ratio between interaction and 
kinetic energies.!^ In ID lattices the HCB Hamiltonian can be mapped onto the ID 
XY model of Lieb, Schulz, and MattisP ^ For perio dic systems, this model has been 
also studied extensively in the literature,l2^EZEHl and more recently in trapsl^As 
in the continuum case, in a periodic lattice the HCB one-particle density matrix 
decays as px ^ l/v^ foi' noninteger fillings, so that the highest occupied effective 
single particle state has an occupation ~ ^/Nt- Remarkably, the above properties 
are universal independently of the power of the confining potential in the trapped 
case, and of the presence of Mott-insulating domains in the system.'^^ 

One fascinating possibility opened by the experimental realization of the TG 
regime is the study of the nonequilibrium dynamics of such strongly correlated 
bosonic systems. As we review in this article, two very peculiar features occur 
during the expansion of a HCB gas on a ID lattice. For low initial densities, the 
momentum dist ribu tion of expanding HCB's rapidly approaches that of noninter- 
acting fermions 11221 This occurs without any evidence of loss of coherence in the 
system, but rather with a slight increment of it, as reflected by the occupation of 
the lowest natural orbital. On the other hand, when the system starts its free evo- 
lution from a pure Mott-insulating (Fock) state, quasicondensates of HCB's emerge 
at finite momentum.!^ 

We also present in this work a detailed study of the out-of-equilibrium behavior 
of the Shannon information entropieJ32in coordinate and momentum space, for the 
two situations mentioned before. As the application of information concepts to the 
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study of correlated systems is increasing, we find that during the expansion of the 
TG gas in a lattice, the Shannon information entropy in momentum space displays 
a qualitative different behavior depending on the initial conditions in the system. 
This behavior allows to identify the regime where a fast fermionization of HCB 
momentum distribution occurs, from the one in which traveling quasicondensates 
are dynamically generated in the system, and study the crossover between them. 

The exposition is organized as follows. In Sec. [21 we briefly describe the exact 
numerical approach we follow to study the nonequilibrium dynamics. We also intro- 
duce concepts that will be used in the rest of the manuscript. In Sec.|31we analyze 
the dynamical fermionization of the HCB momentum distribution function. The 
emergence of quasicondensates at finite momentum is reviewed in Sec. 01 Sectional 
is devoted to the analysis of the Shannon information entropies in coordinate and 
momentum space for the two cases discussed in Secs.Oand01 and their contrast to 
the ones in the fermionic case. Finally, the conclusions are presented in Sec. O 

2. Exact Approach 

The HCB Hamiltonian can be written as 

^ = -^E {blh+,+h.c.)+V2j2^^ '^^ &r = = 0, {b„bl} = 1, (1) 

i i 

where the bosonic creation and annihilation operators at site i are denoted by bl 
and 6j respectively, and the local density operator by = 6|6j. t is the hopping 
parameter, and the last term represents a harmonic trap with curvature V2 and site 
positions Xi — ia {a is the lattice constant). The brackets in Eq. apply only to 
on-site anticommutation relations, for i ^ j these operators commute [b^, b^-] = 0. 

To calculate HCB properties we use the Jordan- Wigner transformation,!^ 

i-1 i-1 

bl = f!Yle~'^fl^f^, b.^l[e'^4Uf, , (2) 

/3=1 /3=1 

which maps the HCB Hamiltonian onto the one of noninteracting spinless fermions 

i i 

where // and fi are the creation and annihilation operators for spinless fermions 
at site i, and n( = f} f.^ is the local particle number operator. 

The above mapping allows to express the equal-time Green's function for HGB's 
in a non-equilibrium system as 

G.,(r) = {^B{rm]\^B{r)) = F{r)\ J] e^^^f, fj]^e-^^f'i-<\^ ^{r)), (4) 

13=1 7=1 

where r is the real time variable, 1^'^ (''")) is the time evolving wave function for 
HGB's and \^'p[t)) is the corresponding one for noninteracting fermions. 
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Since there are no interactions in the equivalent fermionic system, the time 
evolution of an initial wave- function (l^*^)) can be calculated numerically as 

Nf N 

(5=1 CT=1 

which is a product of single particle states, with Nf being the number of fermions 
(Nf = Nb), N the number of lattice sites, and P(t) the matrix of components of 
1^1' i?(T)). Once |^'i?(r)) is calculated, the action of Yij^i'^^^'"'^'^'^'' the right in 
Eq. Q generates only a change of sign on the elements Pas{T) for cr < j — 1, and 
the further creation of a particle at site j implies the addition of one column to 
P(t) with the element PjAr^,_|_i(T) = 1, and all the others equal to zero. The same 

can be done for the action of n^=\ e'''"^^^^'^ fi on the left of Eq. (g)). Hence, the HCB 
Green's function can be obtained aJ 



Nf + l N Nf + 1 N 

Gr,{r) = (oi n E w/^ n E ^7?M/7 io) 

5=1 (3=1 CT=1 7=1 



det 



(P'^(r))'p'^(7 



(6) 



where P' (r) and P' (r) are the new matrices calculated from P(t) when the 
required signs are changed and the new columns added. The last equality in Eq. 
((HJ is obtained after considering the following identity 

(0|/(Ti • • ■ fa-N^ + lf&Nf + l ' ' ■ /ffJO) = £ ^ "^^^^criS-Xi ' ■ ' '^(T]Vj,+l5-Ajy^^j J (7) 

where £^^"'^"f+^ ig the Levi-Civita symbol in Nf + 1 dimensions, the indices A have 
values between one and Nf + 1. 

The evaluation of Gij{T) is done numerically, and the equal-time one-particle 
density matrix [/Oij (T)] is determined by the expression 

P^J{T) = (*B(r)|6,|6,|*B(r)) - G,.(r) + % [1 - 2G.,(r)] . (8) 

There are two quantities related to pij that will be studied here because of 
their relevance to experiments, and to the understanding of the coherence of the 
system. These quantities are the momentum distribution function (rifc), and the 
so-called natural orbitals. Uk is usually obtained experimentally in time-of-flight 
measurements. We calculate it as 

n,{r)^{a/0j2e-"'^'''^PAr), C = • (9) 

Notice that we normalize nk using which is a length scale set by the combination 
lattice-harmonic confining potential in Another quantity that is relevant to 
define by means of ^ is the characteristic density p = It plays for 

trapped systems a role similar to the one of the density p — N^/N \a. the periodic 
case. As p — > 0, the lattice effects dissappear, and one recovers the continuum limit. 



I 



I 



2, 2008 16:2 WSPC/INSTRUCTION FILE Expansionl 



5 

On the other hand, as p increases b eyond a critical value (pc ~ 2.6—2.7) a MI region 
builds up in the middle of the trap.^^HIMl 

Also of interest are the natural orbitals ((/)^). They can be considered as effective 
single-particle states in these strongly interacting systems, and are defined as the 
eigenf unctions of 

N 

having occupations A,j. In dilute higher dimensional gases, when only the lowest 
natural orbital (the highest occupied one) scales ^ iV^, it can be regarded as the 
BEG order parameter, i.e., the condensate.36 Although for periodic systems in 
equilibrium the natural orbitals are momentum states, this is not in general the 
case for trapped or out-of-equilibrium systems. 

3. Dynamical Fermionization 

In Fig. n^a) we show the evolution of Uk for 100 HCB's once the harmonic trap 
confining the system is turned off. At r = 0, the characteristic density in the trap is 
p — 0.51, i.e., far from the regime with a Mott insulator. We compare the HCB Uk 
with the one of the equivalent noninteracting fermions (which remains unchanged 
during the expansion). The Fermi momentum kp is defined as ep = —2tcos{kpa), 
where ep is the energy of the last occupied fermionic single-particle state in the 
trap at r = 0. In contrast to a periodic system, in a trap is continuous aX kp. 
However, it approaches zero even faster than an exponential for k > kp. 




71/8 7t/4 kp37t/8 71/2 1 10 100 1000 

ka ^ 

Fig. 1. nj. (a) and A (b) for 100 HCB's expanding from an initial state with V20? = 2.6 X 10~^t, 
and compared to the ones for the corresponding fermions. Times (r) are given in units of h/t. 
In (a), kp denotes the Fermi momentum, and a the lattice constant. In (b), the thin dashed line 
corresponds to a power law »7~*, which is known from equilibrium systems (see text). 

Two remarkable features are evident in Fig. ^a). (i) Shortly after switching 
off the trapping potential, the peak at nk=Q disappears, (ii) The expansion of the 
system leads to an for the HCB's that is equal to the one of the fermions. 
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This fermionization of Uk starts from the low moraentum states towards kp, and 
produces a Fermi edge. 

After observing the above behavior of the HCB Uk, one could infer that the 
system is loosing coherence, i.e. that the effective single-particle states obtained 
after diagonalizing pij may reduce their occupations (A) with time toward A = 1, 
the value for noninteracting fermions. As shown in Fig. Q^b) this is not the case. 
The occupation of the natural orbitals almost does not change. Actually, as seen in 
Fig. |2Ia) the lowest natural orbital occupation Aq slightly increases its occupation 
instead of reducing it. 

The out-of-equilibrium increase of Aq is similar to one in the ground state when, 
keeping constant the number of particles, the curvature of the trap is decreased.!^ 
In both cases the increment of Ap can be intuitively understood as an enhancement 
of the coherence in the system due to an increase of its size, which delocalize the 
HCB's over more lattice sites. However, in equilibrium, nk=o also increases along 
with Aq. The different behavior of Aq vs nk=o in- and out-of-equilibrium can be 
understood since in the last case the lowest natural orbital is composed by HCB's 
with many different momenta. This can be seen in Fig. ^[h) , where we show the 
Fourier transform of the lowest natural orbital (|</)°|) at different times. At r = 
one can see that |0°| has a peak at fc = showing that quasi-condensation occurs 
around fc = 0, and this is reflected in Uk- For r > the lowest NO Fourier transform 
extends in /c-space so that 0° starts to be composed by HCB's with low and large 
momenta, i.e., its structure changes in momentum space during the expansion. 
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Fig. 2. Lowest natural orbital occupation (a) and the Fourier transform of the lowest natural 
orbital (b) for 100 HCB's expanding from an initial state with V2a^ = 2.6 X lO^^t. Times (r) are 
given in units of h/t. In (b), kp denotes de Fermi momentum. 



In Fig. n^b) there is a further salient feature, which sets in when the density of 
the expanding HCB's becomes very low. A decay A^ ^ ?7~^ develops for large values 
of rj. In equilibrium we have shown that such a power-law deca y is universal at low 
densities, independently of the power of the confining potentiaLl^ The prefactor of 
the power law ^7V(, was found to depend only on A^f,.^^ We find out-of-equilibrium 
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that Anj^ has exactly the same value than in the ground-state!^ For — 100 we 
have plotted A,, = A^^rj^^ in Fig.mb). Presumably, the universal decay A,, ^ 7/"^ 
is related to the singular character of the H CB (5-i nteraction, as the case for the tail 
TT-k \k\~'^ obtained for continuous systemsl ^^ l ^'^ lHowever, analytical proofs of this, 
and the universality of the prefact or - Aat^ , have not been given so far. The power 
law Uk ~ 1^1 that in equilibriurJ^ appears together with A,, ~ ry^'*, disappears 
during the expansion, and the HCB rifc starts to behave like the one of the fermions. 

An overall understanding of the previously discussed nonequilibrium behavior of 
rife, and natural orbital occupations, can be gained directly studying the one-particle 
density matrix. Out of equilibrium p,y = |py |e**»J , i.e., it is complex. Results for 
Pij in the same systems of Figs. ^ and |2Ib) are presented in Fig. 13 Fig.|3{a) shows 
that |pij(T)| exhibits the same power-law decay than in the ground state,!^ 
i.e., \pij\ ~ \xi — Xj\~^/^ for large values of \xi — Xj\ and for all times. Hence, 
this decay of the one-particle correlations is the one accounting for the large, and 
increasing, occupation of the lowest natural orbital as the system expands (like in 
the ground state). On the other hand, Figs.|2tb)-(d) show that the phase of pij {9ij) 
starts to increasingly oscillate at large distances. This phase is the one accounting 
for the fermionization of Uk- In particular. Fig. Ufb) shows that after a very short 
time, when the modulus of the OFDM have almost not changed, 9ij has started to 
oscillate for \xi — Xj \ ^ a producing a fast destruction of the zero momentum peak 
in n^, as shown in Fig.^a). 




\jCi — Xj \/a \Xi — Xj \/a 

Fig. 3. Modulus of pij (a) and its phase (b)-(d) for the same initial trap parameters and times of 
Figs, m and l2lb^. Both quantities have been calculated with respect to the center of the system. 
Thin continuous lines in (a) correspond to power laws \xi — Xj\~^^'^. 



So far we have presented results for a trap with a characteristic density p — 0.51 
and Nh — 100. In what follows we analyze the effects of changing A^f, and p in the 
system. For that, we study the relative area between nfc for HCB's and ri{ for 
fermions, (5 = |rifc — ri{|)/(^^ rifc). We will consider that the Uk of the HCB's 
has fermionized when 8 < 0.05, i.e., when its difference with the one of the fermions 
is less or equal than a 5%. 
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In Fig. Efa) we show (5 vs r in systems with p — 0.51 when the number of 
particles in the trap is increased. Not ice that keeping p constant is equivalent to 
keep constant the Fermi energy e p 1^ Figure Ela) , and its inset, show that in this 
case the fermionization time {tf) increases linearly with Ni,. The fast disappearance 
of the fc = peak in Uk [Fig. ^a)] is reflected in Fig. ^a) by a fast reduction of 5. 
We find that after long times the reduction of 5 is very close to a power law. This 
means that tf depends strongly on the criterion chosen. 

The consequences of increasing p, and hence the Fermi energy, are analyzed in 
Fig.^b). In order to compare systems with different p, i.e., different Uk, we display 
in the inset of Fig. Efb) the ratio R between the size of the cloud when 5 = 0.05, 
and its initial size (before expansion). (R is independent of the number of particles 
for a given value of pP'^) The inset in Fig. Efb) shows that with decreasing p the 
ratio R reduces up to ~ 2.5. For low p, such that the interparticle distance is much 
larger than the lattice spacing, the initial lattice gas is equivalent to the one in 
continuous systems. This means that R 2.5 is relevant when there is no lattice, 
where the fermionization of Uk can be more easily observed. In continuous systems, 
the dynamical fermionization of the HCB nk has been also receiitly studied.l^^EHI 
Interestingly, the fermionization time was found to be Tp ~ l/ei?,'^in contrast to 
our findings in the lattice where it increases linearly with Nh for a fixed ep [inset 
in Fig. El^a)] , and dramatically for large characteristic densities [Fig. 01b)] . 




125 250 375 500 100 200 300 400 



Fig. 4. Fermionization of during the expansion of the gas. (a) Decrease of 5 (see text) as a 
function of time for p = 0.51; A^;, = 51 (dashed line), A^;, = 100 (continuous hne), Ni, = 150 
(dotted hne), and Ni, = 200 (dashed-dotted line). The inset shows the fermionization time Tp (see 
text) vs N),, for p = 0.51. (b) Decrease of 5 as a function of time for Ni, = 100, p = 2.0 (dashed 
line), p = 2.5 (continuous line); and Nj, = 101, p = 3.0 (dotted Une), p = 4.0 (dashed-dotted hne). 
In the last two cases there is a Mott-insulating domain in the center of the trap. We have plotted 
as a thin continuous line 5vs r for Ni, = 101 in a pure Mott insulating state, i.e., a state with one 
particle per lattice site. The inset shows the ratio R between the size of gas when 5 = 0.05 and 
its original size vs p for Ni, = 100. 

As shown in Fig. 2fb), increasing the characteristic density beyond the values 
depicted in its inset, the behavior of 6 starts to depart from the one seen in Fig. 
0Ja). This is because particles become more localized in the middle of the trap. 
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and after p ^ 2.6 — 2.7 a Mott insulator appears in the system. This localization 
generates an Uk that approaches the of the fermions in the initial state. [In the 
limit where all occupied lattice sites have rn — 1, nk{T — Q) — n{.] When such 
systems are released from the trap 5 increases for some time. We will study in the 
next section the origin of this increase. We should mention, however, that after long 
times a fermionization of n^. starts to occur as before for smaller p. The difference 
is that, as shown in Fig.^fb), the time scales for this process are very large. 

4. Emergence of Quasi-Condensates at Finite Momentum 

We study in this section the free expansion of HCB's that are initially in a pure 
Mott-insulating state with one particle per lattice site. Since double occupancy 
is forbidden for very large on-site repulsions U , such states can be created ex- 
perimentally in very deep optical lattices, and using strong confining potentials 
\U » V2(Affca/2)^ » t], so that vacancies in the initial state are suppressed. 

In Fig.Elwe show the evolution of density (a) and momentum (b) profiles of 301 
HCB's after the confining potential is turned off. At r = 0, is flat corresponding 
to a pure Fock state, where particles are localized in real space. The corresponding 
7i{ of the fermions is identical, and does not change with time. Surprisingly, sharp 
peaks appear in the HCB at A; = ±7r/2a during the expansion. Such peaks are 
very similar to the ones present at fc = in the ground state of a system with 
superfluid domains .12^ 




-500 -250 250 500 -tz -7i/2 7t/2 n 



x/a ka 

Fig. 5. Evolution of density (a) and momentum (b) profiles of 301 HCB's in 1000 lattice sites. The 
times are t = (A), 60h/t (0)i 120?i/t (x), and 160/i/t (V). Positions (a) and momenta (b) are 
normalized by the lattice constant a. 

Since the peaks at ni^^j.^/2a indicate the presence of quasicondensates, like 
in the ground state, but at finite momentum in this case, we diagonalize pij and 
study the natural orbitals and their occupations. In Fig.EJa) we show the evolution 
of the lowest 500 natural orbital occupations during the expansion. Like in nk, one 
can see that with the disappearance of the A = 1 plateau, which signals the Mott 
insulating state, a peak appears for the highest occupied orbitals. As seen in the 
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inset of Fig. Efa) (where we only show the eleven highest values of A), the two 
highest occupied orbitals are degenerate. We discuss below that the two peaks seen 
in Uk signal their emergence. 
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Fig. 6. Evolution of the natural orbital occupations (a) and the lowest natural orbital wave function 
(b) of 301 HCB's in 1000 lattice sites. The times are r = (A), 6Qh/t (Q), 120;i/i (x), and 160^* 
(V), and correspond to the density and momentum profiles shown in Fig. 1^ 



The evolution of the modulus of lowest natural orbital wave function {(jp is 
complex) is shown in Fig. IH^b). The results for \(f>^\ are identical to those of |</>°|. 
(At T — the natural orbitals are delta functions at the occupied sites, so that we 
do not show |(/)°(t — 0)| in the figure.) As seen in Fig. |SJb) during the expansion 
of the HCB's, two identical wave packets (or lobes) form at the sides of the Mott 
insulator. Once the Mott-insulating domain has disappeared, the shape of the wave 
packets almost does not change, like for a soliton, and they just move in opposite 
directions. To show that we have superposed in Fig. [Tfa) the right moving lobe at 
three different times. They are almost indistinguishable. This solitonic character 




Fig. 7. (a) Superposed right lobe of the lowest natural orbital at t = 120h/t (x), lAOh/t (A), and 
160/i/t (V). The inset shows the evolution of the lowest natural orbital right lobe's position. The 
line going through the points has a slope 2. (b) Fourier transform of the lowest natural orbital at 
T = 60^* (0)> 120?i/i (x), and 160h/t (V). 



I 



I 



2, 2008 16:2 WSPC/INSTRUCTION FILE Expansionl 



11 

of the moving wave packets allows us to calculate their group velocity, which we 
find to be constant and equal to vno — i'iat/h. This is shown in the inset of Fig. 
Ufa) where we have plotted the position of the right moving wave packet xno vs 
T, in units of h/t, (the results for the left moving one are identical with x — > —x), 
which is a straight line with slope 2. Given the HCB dispersion relation in a lattice 
e/j = —2tcoska (identical to the one of noninteracting fermions) one can obtain, 
using the well known expression v = {l/h){dek/dk), the momenta corresponding 
to the moving lobes to be fc = ±7r/2a. 

The above result can be further confirmed by calculating the Fourier transform 
of the lowest natural orbital at different times. They are depicted in Fig.lJJb), and 
exhibit sharp peaks centered at quasi- momenta k = ±7r/2a. The peak at fc = +7r/2a 
appears due to the Fourier transform of the right lobe in Fig. El^a) (i.e. 0" ~ 
|^0|gi7rxi/2a ^ q-j^ ^^iq quc at fc = —TT /2a duc to the fourier transform of 

the left lobe (i.e. 0^ ~ l^^je"""^'/^" for Xi < 0). We have also calculated the Fourier 
transform of the other natural orbitals finding that they have a very small or zero 
weight at fc = ±7r/2a. Hence, the peaks appearing at fT.fc=±T/2a are reflecting the 
existence of the two traveling, with fc = ±7r/2a, highly occupied natural orbitals. 

Apart from their large occupations, the two lowest natural orbitals exhibit an- 
other remarkable property. Once formed, their moving lobes have exactly the same 
form independent of the number of particles (iVb) in the initial Mott insulating 
state. Their size L is proportional to A^f,. Actually, as shown in Fig.lSfa) the lobes 
can be perfectly rescaled considering \ipo\ ~ A^^^^^j^^l vs x/N^a. 



1 I ■ 1 1 1 1 1 1 F I I b 




x/Ni,a \x,-Xj\/a 

Fig. 8. (a) Scaled right lobe of the lowest natural orbital (see text) vs x/Ni,a for Ni, = 101 (0)i 
201 (x), and 301 (V). (b) One-particle density matrix measured with Xi fbced at the left extreme 
of the left lobe, and Xj > Xi. Times are the same of Figs. ISI7I r=60?!,/t (0)i 120?!,/t (x), and 
IQOh/t (V). The straight line is 0.25 \{xi - Xj)/a\-'^/^ (see text). 

In order to understand the origin of the large occupation of the two lowest 
natural orbitals, and their scaling properties, we study the one-particle density 
matrix. At r = 0, in the Mott insulating state, pij — Sij at the occupied lattice sites, 
and zero in the rest of the system. During the expansion of the HCB's, a surprising 
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behavior sets in. Quasi-long-range one-particle correlations develop in the region 
of the system where the lobes of the lowest two natural orbitals were observed in 
Figs.|nib)-|HIa). The decay of these correlations is \pij\ = 0.25 \{xi — Xj)/a\~^/'^, as 
depicted in Fig. IS^b). They exhibit the same power-law decay that was shown to 
be universal, independent of the exponent of the trapping potential, in the ground 
state.Ei Outside the regions where the lobes are present, pij decays faster as seen 
in the tails of \pij \ present in Fig. IHl^b). These quasi-long-range correlations in two 
disconnected segments of the system are the reason for the degeneracy found in the 
lowest natural orbitals. A detailed analysis of pij also shows that, as expected, the 
peak in Uk &t k ~ +7r/2a is originated by components of pij with Xi,Xj > 0, and 
the one at fc = —t: /2a by the components with Xi, Xj < 0, so that in the regions of 
the lobes one can write pij ~ 0.25 \{xi — a;j)/a|~^/^e*^^^'~^J^/^'^ for Xi,Xj > 0, and 
p,j ~ 0.25 \ix, ~ a;j)/a|-i/2e-"(^'-^^)/2a ^^^^^ ^ q rpj^^ prefactor of the power 

law (0.25) was found to be independent of time, and of the number of particles in 
the initial Fock state. 

Considering all the properties discussed so far one can show that the two high- 
est occupied natural orbitals, dynamically generated and traveling with momenta 
ib7r/ 2 a re quasicondensates in the same sense than those observed in the ground 
state.l^ They have a diverging occupation Aq as the number of particles in the 
initial Mott insulating state increases, however Xq/Ni, as Ni, ^ oo. Aq can be 



obtained as 



Ao = ~ £ dx J^"^ dy 

Jrf.x/'.. l^°'gl°-y/'-)l ,.V^. (U) 

The sums were replaced by integrals considering that L 3> a. In the second line of 
Eq. Hll|l we changed variables x=XNi,a, y=YNi,a, (jP=Nj^ ^^'^ip'^ , and noticed that 
the phase factors between the natural orbital and pij cancel out. The integral over 
X,Y \s & constant that we call A. The final result in Eq. pi|l shows that, like in the 
ground state, Aq ~ ^/T%- This result was confirmed by our numerical calculations 
as shown in Fig. Efa) . There we have plotted the maximum occupation of the lowest 
natural orbital vs N},. A fitting of the constant A reveals a value of 0.72. 

The appearance of two traveling quasicondensates, at momenta k — ±7r/2a, can 
be understood on the basis of total energy conservation. Given the HCB dispersion 
relation in a lattice tk = — 2tcosfca, we notice that the total energy of the initial 
Mott insulating state is zero [Et = 0), since nk is totally flat. Considering that the 
energy is conserved in the system, if all the particles would condense into one state, 
it would be the one with an energy corresponding to Zk — Et/N = 0. For the case 
of free expanding HCB's there are two states with Zk = 0. They are the momentum 
states with fc = ±n/2a. Actually, in ID there is only quasicondensation, so that 
the argument above applies only to maximize the occupation of the fc — ±7r/2a 
states. Two other important properties of the k = ±7r/2a wave packets are that 
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they travel with the maximum velocity possible in a lattice, and that the density 
of states of the free system has minima at these values of k. The latter property 
reinforces quasicondensation into the two single k = ±7r/2a states. 




1/2 

Fig. 9. (a) Maximum occupation of the lowest natural orbital vs A^;,. The straight line is 0.72 N^' . 
The inset shows the time at which the maximum occupation of Aq is reached Tm (in units of h/t) 
vs Nf,. The line following the points is 0.32Afi,. (b) Time evolution of the lowest NO occupation 
for Nt = 101 (O), 201 (x), and 301 (V). The straight line is 1.38y^Tt/h. 

During the formation of the quasicondensates, the increase of their occupation 
is also characterized by a power law. This is shown in Fig.|UJb), where we plot Ao as 
a function of the evolution time. The log-log scale evidences the fast increase of Xq. 
A detailed examination shows that the population of the quasicondensate increases 
in a universal way [LSSy/rt/h, continuous line in Fig. E^b)] independently of the 
initial number of particles in the Fock state. This power law is determined by the 
universal behavior of the off-diagonal part of ptj shown before, and by the fact 
that during the formation of the quasicondensates, they increase their sizes linearly 
with time at a rate given by the maximum velocity in the lattice Ii^atoI = 2at/h. 
The power law mentioned before is followed up to the point where the maximal 
occupation is reached, which corresponds to the point at which the initial Mott 
insulator has completely melted. One can then realize that given vno, the time 
at which this occurs depends linearly on the number of particles. We find that 
Tra — 0.32Ni,h/t, as shown in the inset of Fig.Efa). To have an example of the time 
scales we are referring here, we consider a tipical experimental setup of rubidium 
atoms, in a lattice with a recoil energy of Er = 20kHz, and a depth of 2QEr. In 
these systems, the time Tm can be estimated as r™ ~ 5.7iVb(ms). 

A second characteristic of these systems that can be seen in the log-log plot 
in Fig. inib) is the slow reduction of the lowest natural orbital occupation, with 
the consequent reduction of the nf^^^^^2 peaks. This slow reduction is consistent 
with the slow fermionization of Uk seen in Fig. ^b), and suggest that such an 
arrangement could be used to create atom lasers with a wavelength that can be fully 
controlled given the lattice parameter a. No additional effort is needed to separate 
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the quasi-coherent part from the rest since the quasicondensate is travehng at the 
maximum velocity on the lattice so that the fr'ont part of the expanding cloud is 
the quasi-coherent part. 

The actual experimental realization would imply to restrict the evolution of 
the initial Mott state to one direction, as shown in Fig. EH There we display the 
density and momentmn profiles of 150 HCB's restricted to evolve to the right in 
1000 lattice sites at the same evolution times of Fig. El This figure shows that the 
values of ni^^T^/2{T) are almost the same than the ones in Fig. [3 although the initial 
Fock state has half the number of particles. This is because in this case only one 
quasicondensate is created, so that the lowest natural orbital also has the same 
occupation than in Fig. Efa), and is not degenerate anymore. 




250 500 750 1000 -n -nil 7t/2 u 



x/a ka 

Fig. 10. Evolution of density (a) and momentum (b) profiles of 150 HCB's evolving only to the 
right (the system is open in the left extreme) in 1000 lattice sites. The times are r = (A), 60h/t 
(O), I20h/t (x), and WOH/t (V). 



The above results suggest how to proceed to obtain lasers in higher dimensional 
systems where real condensation can occur.EI One could employ Mott-insulating 
states with one particle per lattice site created by a very strong on-site repulsive 
potential U. The latter is required in order to obtain a close realization of a pure 
Fock state, since quantum fluctuations of the part icle numb er present in a Mott 
insulator for any finite U will be strongly suppressedE32111l Then the geometry of 
the lattice should be designed in order to restrict the evolution of the Mott insulator 
to one direction only, and to have a low density of states around the mean value 
of energy per particle in the initial state. With these conditions the sharp features 
observed in ID may be reproduced by a condensate in higher dimensions. 



5. Information Entropies 

We study in this section the Shannon information entropieJ^ in coordinate (Sn) 
and momentum (Sk) space for the two physical situations discussed in Sees. |21 and 
0] and in between. 
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Here we will follow the definitions 

S'„ = -^riilnni, (12) 

i 

and 

S'fc = - ^ rife In rife, (13) 

k 

normalizing Hi and Hk such that X^j"'* — X^i"'* = 1- Since we have studied the 
expansion of HCB's in systems with ^ 2000 lattice sites, in all the cases discussed 
in this section we discretize the Brillouin zone with 2000 points to calculate Sk- 
in continuous systems, the sum of Sn and Sk satisfies several interesting prop- 
erties. It was rigorously proven by Bialynicki-Birula and Mycielskj^S. that when 
density and momentum distributions are normalized to unity, 

5T-5„ + 5fe>d(l + ln^), (14) 

with d the space dimension. At the Hartree-Fock level, it has been shown that in 
the ground state of atoms St acquires its minimum value.ESI And more recently, 
it has been proposed that 5*7- can be used as a measure of correlation in atomic 
systems 

In what follows, the comparison of Sn and Sk with their fermionic counterparts 
defined in terms of the fermionic densities (5*^ = — Xi^^/^^'^f) ^"^^ momentum 
distributions {S^ = —^k^k^^^k)^ ^^^^ shown to be very useful. Since density 
profiles are identical for HCB's and noninteracting fermions, at any time, S'„(r) 
and SI{t) are the same. This identity does not hold in momentum space. In the 
ground state, given the large occupation of the low momentum states in the TG 
gas, in contrast to the extended rik of the fermions, one can realize that SI{t = 
0) > Sk{T — 0), where the equality only holds in the limit in which HCB's and 
fermions are localized in real space (the systems in Sec.0}. 

We first study the behavior of Sn and Sk for the cases in which a fast fermion- 
ization of Uk was observed in Sec. 01 In Fig. IllT a) we show the evolution of Sn for 
systems with the same characteristic density and different Ni,, like in Fig.^fa). As 
the number of particles in the trap is increased keeping p constant, the initial sys- 
tem size increments proportionally to iV;,, producing an increase of Sn- In order to 
compare 5„ for different number of HCB's we have displaced them to coincide with 
Sn for Nb — 100. This can be easily achieved considering S'n = Sn + ln(100/A^b). 
Fig. Illf a) shows that rescaling t by the number of particles, the time evolution 
of the different systems collapse over the same curve which, as expected, shows an 
increase of Sn as the system expands. 

The behavior of Sk, for the systems in Fig. Ilir a). is shown in Fig. Illf b'). Since 
the characteristic density is kept constant, at t — 0, Sk for different Ni, has the 
same value. During the expansion one can see that Sk overcomes S^ (that does 
not change since the fermions are noninteracting). This means that, taking into 
account that the Shannon entropies in coordinate space are identical for HCB's 
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and fermions, as the TG gas expands St of the HCB's becomes larger than the 
corresponding one 5^ of the fermions. For large times St —> S^ from above. As 
discussed below, the increase of St beyond the value of S*^ is a characteristic of the 
fermionization process studied in Sec. 13 




Fig. 11. Evolution of S'„ (a) and Sk (b) in systems with 2000 lattice sites. The initial state has 
p = 1.0 in all cases, and A^;, = 51 (dashed line), Ni, = 100 (continuous line), Nf, = 150 (dotted 
line), and A^;, = 200 (dashed-dotted line). In (a) we have displaced the Sn of Ni, = 51, Ni, = 150, 
and Ni, = 200 by a constant (see text), to superpose it with the one of A^;, = 100. The insets in (a) 
and (b) show the sum and difference, respectively, of Sn and Sk for Ni, = 100. The thin dashed 
line in (b) shows S^, which does not change in time. 



We have also plotted in the insets of Figs. llir a'l and (b) the sum and differences, 
respectively, of Sn and Sk- The sum shows that St has its minimum value in the 
ground state. Hence, the expansion only produces its increment. The difference 
Sn — Sk in the inset of Figs. llir b') reflects the fast destruction of the nk=o peak at 
very short times, which produces changes in 5*^ that are larger than the ones in Sn, 
so that for short times Sk has the largest contribution to the increase of St- For 
large expansion times Sn is the dominant contribution to St as 5*^ Sj, slowly. 

The properties of the information entropies that we have discussed so far are 
generic for all systems that start their evolution from a state with low density in the 
lattice, and in which a fast fermionization of Uk occurs. Hence, we infer that they 
will also be valid in the continuum casc,^ '^^ which can be obtained as the low 
density limit in the lattice. As the characteristic density in the system increases, 
and a Mott insulator appears in the middle of the trap, the behavior of Sk and the 
way in which Sk Sl change qualitatively. 

In Fig. El we show the time evolution of 5„ and Sk for a system with p = 3.0, 
which at r = has a Mott insulating domain in the middle of the trap [Fig. llbf a')]. 
During the expansion of the HCB's for this initial condition 5„ [Fig. ^Ja)] is very 
similar to the one in Fig. lllf a) for low initial densities, i.e., no distinguishing feature 
appears in Sn for different initial densities. 

On the contrary, the behavior of Sk in Fig. Il^r b) is very different to the one seen 
in Fig. llir b). with the exception of very short times where Sk in Fig. ll^r b) exhibits 
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a small increase. It corresponds to the fast reduction of the nk=a peak similarly to 
the one seen in Fig. ^a). At r = 0, there is a peak a,t k — [Fig. I13f b)] because 
of the superfluid domains surrounding the Mott insulator in the middle of the trap 
[left inset in Fig. llcSr a)] . At short times these superfluid regions are affected by the 
expansion like the systems analyzed in SecOl producing an increase of Sk- 




With the melting of the Mott insulator in the middle of the trap, quasi-coherence 
builds between initially uncorrelated particles producing two peaks in n^. As shown 
in Fig. Il^r bl for p = 3, they are even larger than the one in the ground state. This 
causes the reduction of Sk seen in Fig. I12r b'). Notice that the situation is similar 
to the one discussed in Sec. ^ but now the initial condition is different since the 
Mott insulator is surrounded by superfluid phases. Because of that, the peaks in 
Uk do not appear at fc = ±7r/2 [Fig.^Jb)], although with time they move toward 
k = ±7r/2. In this case the lowest natural orbitals are, during the expansion, more 
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Fig. 13. Evolution of density (a) and momentum (b) profiles of 101 HCB's in 2000 lattice sites, 
for and initial p = 3.0. The times are t = (continuous line), Sih/t (dotted line), 400/i/t (dashed 
line). The density profiles at r = and r = 3Ah/t can be better seen in the insets in (a). In (b) 
we have plotted as a thin continuous line the momentum distribution function of the fermions. 
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complicated than the ones described in Sec. 01 since HCB's with many different 
momenta are couple into them. Hence, in experiments, the quality of the k = ±tt/2 
traveling quasicondensates needs to be ensured by reducing the size of the superfluid 
domains surrounding the Mott-insulating core. 

After the peaks in Uk have acquired their maximum value, Sk starts to slowly 
increase towards n^. As shown in the inset of Fig. I12r a) the sum of 5„ and Sk 
always increases during the expansion. Even when Sk decreases, the increase of 
Sn is larger and warrants that the minimum of St occurs in the ground state. As 
signaled by arrows in Fig. 1121 a kink can be observed in St when 5*^ attains its 
minimum value. 

We conclude this section with the analysis of the information entropies when the 
system starts its evolution from pure Mott-insulating states, like the ones discussed 
in Sec. 01 As shown in Figs. I14f a').fb). and the inset in (a), the behavior of Sn 
and Sk is the same observed in Fig. 1121 [if one excludes the very short times when 
Sk increases in Fig. I12r b) due to the superfluid domains that are not present in 
the case presented in Fig. . Since 5*^ in these systems is always smaller than 
Sl, we have found a criterion in terms of the Shannon information entropy in 
momentum space that allows to distinguish between (i) systems in which a fast 
fermionization of occurs (Sec.|31) where 5*^ S^ from above, and (ii) systems in 
which quasi-coherence develops dynamically producing traveling quasicondensates 
(Sec. 0} where Sk — > Sl from below. The latter feature is more difficult to see for 
the number of particles we have employed all over this work {Nt, = 100-300), due to 
the large system sizes one requires to calculate. However, in the inset in Fig. I14f b) 
we present results obtained for the expansion of 11 HCB's in up to 3000 lattice 
sites were one can clearly see that Sk is always smaller than Sj^ and that Sk —> S[ 
for very long times. 




100 200 300 400 100 200 300 400 



Fig. 14. Evolution of Sn (a) and Sk (b) in a system with 101 HCB's, 2000 lattice sites, and initially 
in a pure Mott insulating state. In the inset in (a) we show the sum of Sn and S^- In the inset in 
(b) we have plotted Sk vs t for an initial Mott insulating state with 11 HCB's and expanding in 
a system with up to 3000 lattice sites. The thin dashed lines show S^ . 



2, 2008 16:2 WSPC/INSTRUCTION FILE Expansionl 



19 

6. Conclusions 

In this work we have presented a unified study of the expansion of TG gases in ID 
optical lattices.'SQ'ffl We have employed an exact numerical approach that allows 
to study systems with large number of particles and sizes, and for arbitrary long 
time scales. We find that two very different regimes can be observed according to 
the initial conditions of the trapped gas. 

(Regime I) When the TG gas starts it free expansion from a low density sys- 
tem, the momentum distribution of the bosons rapidly approaches the one of the 
equivalent noninteracting fermions. This fermionization of Uk occurs without loss 
of coherence in the system. Actually, coherence increases as shown by the increase 
of the occupation of the lowest natural orbitals. This can be understood due to 
the presence of quasi-long range one-particle correlations, which have the same 
power-law decay as in the equilibrium case. A new feature that appears during the 
expansion of the gas is that the lowest natural orbital starts to be populated by 
particles with many different momenta, in contrast to the ground state where it is 
mainly populated by particles with fc = 0. 

(Regime II) When the TG gas starts it free expansion from a high density 
(Mott-insulating) trapped system, quasi-long range correlations develop between 
initially uncorrelated particles. They exhibit the same power-law decay known from 
the ground state, and produce the emergence of quasicondensates of HCB's at fi- 
nite momentum. Their maximum occupation was found to be proportional to the 
squared root of the number of particles in the system. This quasi-condensation 
out-of-equilibrium is reflected by the appearance of two identical peaks in the mo- 
mentum distribution function. We find that the dynamically generated quasicon- 
densates are very similar to solitons since they can travel long distances almost 
without changing the shape of the wave-function and their occupation. Their mo- 
menta can be fully controlled by means of the lattice parameter. 

Finally, we have shown that the Shannon information entropy in momentum 
space provides a criterion to distinguish the two regimes previously mentioned. In 
the first regime the Shannon information entropy in momentum space becomes 
larger than the corresponding one of the fermions, and with time approaches the 
last one from above. For the second regime we find that the Shannon informa- 
tion entropy in momentum space is always smaller for HCB's than for fermions so 
that the first one approaches the second one from below. (The fermionic entropy 
in momentum space remains unchanged during the free expansion since they are 
noninteracting particles.) The information entropy in momentum space was also 
shown to be useful to study the crossover between both regimes. It allows to dis- 
tinguish the times scales in which the different processes related to the superfluid 
and Mott-insulating domains set in the momentum distribution function. 



February 2, 2008 16:2 WSPC/INSTRUCTION FILE Expansionl 



20 

7. Acknowledgments 

We are grateful to R. T. Scalettar and R. R. P. Singh for stimulating discussions. 
This work was supported by NSF-DMR-0312261, NSF-DMR-0240918, NSF-ITR- 
0313390, and SFB 382. We thank HLR-Stuttgart (Project DynMet) for allocation 
of computer time. 

References 

1. F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, T. Bourdel, J. Cubizolles, and C. 
Salomon, Phys. Rev. Lett. 87 (2001) 080403. 

2. A. Gorlitz, J. M. Vogels, A. E. Leanhardt, C. Raman, T. L. Gustavson, J. R. Abo- 
Shaeer, A. P. Chikkatur, S. Gupta, S. Inouye, T. Rosenband, and W. Ketterle, Phys. 
Rev. Lett. 87 (2001) 130402. 

3. T. P. Meyrath, F. Schreck, J. L. Hanssen, C.-S. Chuu, and M. G. Raizen, Phys. Rev. 
A71 (2005) 041604(R). 

4. M. Greiner, I. Bloch, O. Mandel, T. W. Hansch, and T. Esslingcr, Phys. Rev. Lett. 87 
(2001) 160405. 

5. H. Moritz, T. Stoferle, M. Kohl, and T. Esslinger, Phys. Rev. Lett. 91 (2003) 250402. 

6. T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. Esslingcr, Phys. Rev. Lett. 92 (2004) 
130403. 

7. B. L. Toha, K. M. O'Hara, J. H. Huckans, W. D. Phillips, S. L. Rolston, and J. V. 
Porto, Phys. Rev. Lett. 92 (2004) 190401. 

8. B. Paredes, A. Widera, V. Murg, O. Mandel, S. FoUing, I. Cirac, G. V. Shlyapnikov, 
T. W. Hansch, and I. Bloch, Nature 429 (2004) 277. 

9. T. Kinoshita, T. Wenger, and D. S. Weiss, Science 305 (2004) 1125. 

10. C. D. Fertig, K. M. O'Hara, J. H. Huckans, S. L. Rolston, W. D. Phillips, and J. V. 
Porto, Phys. Rev. Lett. 94 (2005) 120403. 

11. M. Girardeau, J. Math. Phys. 1 (1960) 516. 

12. M. Olshanii, Phys. Rev. Lett. 81 (1998) 938. 

13. D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 85 (2000) 
3745. 

14. V. Dunjko, V. Lorent, and M. Olshanii, Phys. Rev. Lett. 86 (2001) 5413. 

15. E. H. Liob and W. Liniger, Phys. Rev. 130 (1963) 1605. 

16. A. Lenard, J. Math. Phys. 5 (1964) 930; 7 (1966) 1268. 

17. H. G. Vaidya and C. A. Tracy Phys. Rev. Lett. 42 (1979) 3. 

18. F. D. M. Haldane, Phys. Rev. Lett. 47 (1981) 1840. 

19. V. E. Korepin, N. M. Bogoliubov, and A. G. Izergin, Quantum inverse scattering 
method and correlation functions (Cambridge University Press, Cambridge, New York, 
1993). 

20. M. A. Cazalilla, J. Phys. B37 (2004) SI. 

21. M. D. Girardeau, E. M. Wright, and J. M. Triscari, Phys. Rev. A63 (2001) 033601. 

22. T. Papenbrock, Phys. Rev. A67 (2003) 041601(R). 

23. P. J. Forrester, N. E. Frankel, T. M. Garoni, and N. S. Witte, Phys. Rev. A67 (2003) 
043607. 

24. D. M. Gangardt, J. Phys. A37 (2004) 9335. 

25. E. Lieb, T. Shultz, and D. Mattis, Ann. Phys. (NY) 16 (1961) 407. 

26. B. M. McCoy, Phys. Rev. 173 (1968) 531. 

27. H. G. Vaidya and C. A. Tracy Phys. Lett. A68 (1978) 378. 

28. B. M. McCoy, J. H. Perk, and R. E. Schrock, Nucl. Phys. B220 (1983) 35; B220 



2, 2008 16:2 WSPC/INSTRUCTION FILE Expansionl 



21 

(1983) 269. 

29. M. Rigol and A. Muramatsu, Phys. Rev. A70 (2004) 031603(R); A72 (2005) 013604. 

30. M. Rigol and A. Muramatsu, Phys. Rev. Lett. 94 (2005) 240403. 

31. M. Rigol and A. Muramatsu, Phys. Rev. Lett. 93 (2004) 230404. 

32. C. E. Shannon, Bell Syst. Tech. J., 27, 379 (1948); 27 623 (1948). 

33. P. Jordan and E. Wigner, Z. Phys. 47 (1928) 631. 

34. M. Rigol and A. Muramatsu, Phys. Rev. A7G (2004) 043627. 

35. O. Penrose and L. Onsager, Phys. Rev. 104 (1956) 576. 

36. A. J. Leggett, Rev. Mod. Phys. 73 (2001) 307. 

37. A. Minguzzi, P. Vignolo, and M. P. Tosi, Phys. Lett. A294 (2002) 222. 

38. M. Olshanii and V. Dunjko, Phys. Rev. Lett. 91 (2003) 090401. 

39. B. Sutherland, Phys. Rev. Lett. 80 (1998) 3678. 

40. A. Minguzzi and D. M. Gangardt, Phys. Rev. Lett. 94 (2005) 240404. 

41. E. H. Lieb and R. Seiringer, Phys. Rev. Lett. 88 (2002) 170409. 

42. G. G. Batrouni, V. Rousseau, R. T. Scalottar, M. Rigol, A. Muramatsu, P. J. H. 
Denteneer, and M. Troyer, Phys. Rev. Lett. 89 (2002) 117203. 

43. M. Rigol, A. Muramatsu, G. G. Batrouni, and R. T. Scalettar, Phys. Rev. Lett. 91 

(2003) 130403. 

44. M. Rigol and A. Muramatsu, Phys. Rev. A69 (2004) 053612; Opt. Commun. 243 

(2004) 33. 

45. I. Bialynicki-Birula and J. Mycielski, Commun. Math. Phys. 44 (1979) 129. 

46. S. R. Gadre, S. B. Sears, S. J. Chakravorty, and R. D. Bendale, Phys. Rev. A32 
(1985) 2602. 

47. N. L. Guevara, R. P. Sagar, and R. O. Esquivel, Phys. Rev. A67 (2003) 012507. 



